LOCAL UNIVERSALITY OF REPULSIVE PARTICLE SYSTEMS 
AND RANDOM MATRICES 



FRIEDRICH GOTZE AND MARTIN VENKER 

Abstract. We study local correlations of certain interacting particle systems on the real 
line which show repulsion similar to eigenvalues of random Hermitian matrices. Although the 
new particle system does not seem to have a natural spectral or determinantal representation, 
the local correlations in the bulk coincide in the limit of infinitely many particles with those 
known from random Hermitian matrices, in particular they can be expressed as determinants 
of the so-called sine kernel. These results may provide an explanation for the appearance 
of sine kernel correlation statistics in a number of situations which do not have an obvious 
interpretation in terms of random matrices. 



1. Introduction and Main Results 

This paper is motivated by the surprising emergence of sine kernel statistics in many real 
world observations, parking cars, perching birds on lines and so on. In the field of random 
matrices, the sine kernel describes the local correlations of eigenvalues in the bulk of the 
spectrum of Hermitian random matrices. There it has been shown to be universal to a high 
extent, i.e. it appears for many essentially different matrix distributions. In this article we 
show that the sine kernel describes the local correlations of more general repulsive particle 
systems on the real line which only share the repulsion strength exponent (3 = 2 with the 
eigenvalues of (unitary invariant) Hermitian random matrices. We expect that this behavior 
extends to larger classes of invariant ensembles of random matrices, with repulsion exponents 
(3 different from two. 

To formulate our results, let us recall the so-called invariant /3— ensembles from random matrix 
theory. Given a continuous function Q : M — > R of sufficient growth at infinity and /? > 0, 
set 

PN,QA-)--=^T{\-^--^'e-^^UQ^^^\ (1) 

'^N,Q,I3 J- 

(With a slight abuse of notation, we will not distinguish between a measure and its density.) 
For the "classical values" /3 = 1,2,4, Pn,q,i3 is the eigenvalue distribution of a probability 
ensemble on the space of {N x A^) matrices with real symmetric (/3 = 1), complex Hermitian 
{(3 = 2) or quaternionic self-dual (/3 = 4) entries, respectively. For arbitrary /3, only for 
quadratic Q, Pn,q,(3 is known to be an eigenvalue distribution. 

The notion of bulk universality is usually formulated via the correlation functions of the 
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ensemble. For a probability measure P]\f{x)dx on ! 
function p'^ : R'^ — > M of P/v is defined as 



and A; = 1, 2, 



, N, the k-th correlation 



Pjv(x)(ixjt_|_i . . . dxjv. 



The correlation functions are the densities of the marginals of P/y. The measure p'^{t)dt 
on M.^ is called fc-i/i correlation measure. 

It is known that under very mild conditions on Q, there is an absolutely continuous probability 
measure pQ^fs{t)dt on M, which is the weak limit of p\f q i3{t)dt as N ^ oo. 
Now, Pn,q,i3 is said to admit bulk universality, if for all a with pQ^p[a) > and all ti, . . . ,tk 
the limit 

lim s. p%(a + , ,a+ , , ) (2) 

N^^ PQ^pia)'''^'' \ NpQ^^iaY NpQ^p{a)) 

exists and coincides with the one for Pn,g,I3-, G quadratic (the so-called Gaussian /3-ensemble) . 
Universality here should be understood as coincidence of the limit ([2]) with the corresponding 
Gaussian /3-ensemble. This has been established for large classes of Q. The scaling in ([2]) is 
chosen such that the asymptotic mean spacing between consecutive eigenvalues is normalized 
to 1. However, it is known that the limit depends on /? though. 

In the case /3 = 2, which appears frequently in "real world statistical studies", the limiting 
object ^ is determinantal of type 



lim 



\kPN 



a + 



involving the sine-kernel 



,a + 



^/^Q,2(a) 



det 



sin {'K{ti - tj)) 
■n{ti - tj) 



(3) 



l<i,j<k 



, sinfvrt) , , 
§(t) :=^^,f^0, S(0) :=1. 

TTt 

Universality for unitary invariant ensembles, i.e. (3 = 2 invariant ensembles, was proved in 
many papers, for example (naming only few) [PS97l iPSOSl lDKM+99[ rMMOSi ILLOS] . Recently 
universality (for general /3— Ensembles) was proved in |BEY1H IBEY12] . For /3 = 1,2, bulk 
universality was also proved for Wigner matrices by two groups of authors. Based on earlier 
work of Johansson |Joh01| , universality was shown for general classes of Wigner matrices in 
a series of papers by Erdos, Yau, Schlein, Yin, Ramirez and Peche (see [EY12| for a survey 
on their results) and Tao and Vu (see [TV 12] for a survey on their results). We remark that 
bulk universality was proved in |GG08| for the Hermitian fixed trace ensemble, a random 
matrix which is neither a Wigner matrix nor unitary invariant. 
Writing the density ([T|) in the Gibbsian form 

P^^Q^p = _^e/3E,<,iog|x.-x,|-^Ef=iQ(-.)^ (4) 



Zn,Q,I3 

we see that Pn,q,i3 can be interpreted as an interacting particle system on M in an external 
field, interacting via a 2d Coulomb potential. 

It is believed that many complicated, strongly correlated systems share the local bulk scaling 
limit (defined again by correlation functions) with some random matrix model. This was 
conjectured by Wigner who used random matrices to model energy levels of nuclei. By the 
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underlying matrix structure, physical requirements (conserved quantities, time reversal,. . . ) 
determine the value of /3 in the cases /3 = 1, 2, 4. The limits with /3 = 2 also seem to appear 
in statistics of distances between parking cars [AM06j . waiting times at bus stops in certain 
cities [KSOObj (see [BBDS06] for a determinantal model) and the pair correlation conjecture of 
Montgomery |Mon73j for the zeros of the Riemann Zeta function on the critical line. See e.g. 
[KSOOa] for more relations between the Riemann Zeta function and random matrix theory. 
A common cause for the appearance of sine kernel statistics in a number of statistics about 
real world repulsive systems and in physics and mathematics still remains to be identified. 
We consider here a class of more general interacting particle systems, defined by the density 

^l[^ix,-x,)e-'^^U^(^^\ (5) 

where Q is a continuous function of sufficient growth at infinity compared to the continuous 
function ip : M — > [0, oo). Apart from some technical conditions we will assume that 

(n(0) = 0, ip(t) > for t / and lim ^ = c> 0, (6) 
or, in other terms, is the only zero of and it is of order f3. 

We expect that (at least under some smoothness and growth conditions) the bulk scaling limit 
of ([5]) coincides with that of the /3-ensembles, since in view of the regular local distribution 
of eigenvalues/particles at spacings only the exponents of the interaction kernel should 
determine the local universality class. 

The purpose of this paper is to prove this for /3 = 2 and a special class of (p and Q. From now 
on, we will always deal with the case /3 = 2, therefore omitting the subscript f3. To state our 
results, let h he a continuous even function which is bounded below. Let Q be a continuous 
even function of sufficient growth at infinity. By Pj^ q we will denote the probability density 
on defined by 

1 ^ 
Pn,q(^) := ^ n \xi - exp{-N Q{xj) - ^ hixi - x,)}, (7) 

i<j j = l i<j 

where Z'^ q denotes the normalizing constant. The density -PjvQ ^^^^ written in the 
form dS]) with ip{t) := t'^ exjp{—h{t)}. The first result describes the global scaling limit of 
the correlation measures of Pj^ q- To formulate it, introduce for a twice differentiable convex 

function Q the quantity aq := mit^^Q" {t). Moreover, denote by p^fq the A:— th correlation 
function of P^ q. 

Theorem 1. Let h be a real analytic and even Schwartz function. Then there exists a 
constant > such that for all real analytic, strictly convex and even Q with aq > a^, the 
following holds: 

There exists a compactly supported probability measure /ig having a non-zero and continuous 

density on the interior of its support and for k = 1,2, ... , the k-th correlation measure of 

h f h\'^^ 

P^ Q converges weakly to the k-f old product ifJ-Qj , that is for any bounded and continuous 
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function g : 



lim igpf^QdH= gdU) . (8) 



Af-5-oo _ 

Remark 2. 

a) If h is (additionally) positive definite, then in Theorem [T] may be explicitly chosen 
as = supigK —h"{t). 

b) In general, the measure jJg depends on h. 

c) Q does not seem to be either determinantal or have a natural spectral interpretation, 
therefore we will speak of particles instead of eigenvalues. 

The next result states the universality of the sine kernel in the local scaling limit in the bulk. 

Theorem 3. Let h and Q satisfy the assumptions of Theorem[l\ Then for k = 1,2, . . . , we 
have 



y h,k , h , \ , . 

lim , , , , pj.r n a H , , , , .... a H , , , = det 



sin (7r(ti -tj)) 
7r{ti - tj) 



(9) 

l<i,j<k 



uniformly in ti, . . . ,tk from any compact subset ofK^ and uniformly in the point a from any 
compact proper subset of the support of fiQ. 

Remark 4. If h is positive definite, then in Theorem [3] may explicitly chosen as = 
suptei;-/i"(t). 

We shall demonstrate our approach to bulk universality by means of the following very 
example of functions h and Q. 

Theorem 5. Letj > and a > be arbitrary. Let h{t — s) := 7(t — s)^ andQ{t) = at^ . Then 
([8]) and ([9]) hold for {P^Q^^ uniformly as in Theorem\^ Here Hq will be the semi-circle 
distribution with support [—u},u}], u := (-y/a~+^) ^. 

A first step in the proof of Theorems [1] and [3] is to compare the correlation functions of q 
with correlation functions of eigenvalues of some unitary invariant ensemble. To construct 
such an ensemble, we first determine as the equilibrium measure of some external field V 
(depending on h and Q) using a fixed point argument. The difference between -P^q and this 
unitary invariant ensemble Pn,v consists of (up to normalization) a factor exp{U{x)}, where 
U is a quadratic interaction energy which may be expressed as a mixture of linear interaction 
energy terms using Gaussian processes. This finally leads after some truncation procedure to 
a mixture representation of P^ q by invariant ensembles with the same bulk universality. 
The paper is organised as follows. In Section 2, the asymptotics of Pjy q for h{t — s) := 
'y{t — s)^ and Q{t) = at^ are investigated, in particular Theorem [5] is proved. In Section 
3, we associate to -PJvq a unitary invariant ensemble which will turn out to have the same 
asymptotic behaviour as Pjq q- Section 4 contains concentration of measure inequalities 
which are needed furtheron. Section 5 deals with bounds on the first correlation function 
of a unitary invariant ensemble. The proofs in this section use established techniques which 
we decided to include in detail for the sake of completeness of the exposition. Theorems [1] 



and [3] are proved in Section 6. In the appendix we recall a number of results on equilibrium 
measures. 

A prior version of these results is based on the PhD-Thesis of the second author [Venll] . 

2. A FIRST Example 
In this section, we will study the probability measure 

Pn^{x):=^^ \xi- Xj\^ eyi^{- aNM2{x) --f^{xi- Xjf], (10) 

^ ^<i<j<N i<j 

using the potentials Mp{x) := Xlj^i with p = 2 and constants a, 7 > 0, where Z'^''^ denotes 
the normalization factor. In the following we shall suppress the dependencies on a and 7. 
We will reduce bulk universality of {P^'^) to the well-known bulk universality of the GUE. 
It is convenient to introduce the distribution GUE^j, depending on a parameter a; > 0, as 

PNl^'i^) ■■= i n kfc -^il'expj - AiVM2(x)}. 

Under this scaling the first correlation measure of P^^F will converge to the semicircle law 
supported on [— (for a proof see e.g. |Pas99) ). First we rewrite the density P^r := P^'"' 
using 

7^{xi - Xjf = -fNM2{x) - -fMi{xf, as 

Pn{x) = ^ W -x,f exp{ -(a + 7)7VM2(x)+7Mi(x)2}. (11) 
Using the simple identity 

exp{7t^} = 7;- I exp {e-^t} exp { — e^/4}(ie, we may write (12) 
1 f 1 

^iv(x) = — / — Yl -Xjf exp{ - (a + 7)iVM2(x) + ^eAfi(a;)}exp{-eV4}c^e, 
= ^ [ ^PUx)e-'''^de, where (13) 



Pn{x):=-^ n ki-Xjf exp{-(a + 7)A^M2(x) + ^eMi(x)}, 

■= IT \xi-Xj\^exp{-{a + 'y)NM2ix) + ^eMi{x)}dx. 



l<i<j<N 



We have thus expressed P/v as a probabilistic mixture of the probability measures Pf 



The next lemma deals with the ratio Zpf/Z^. 
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Lemma 6. For each e, each N and all 0,7 > we have 



.2 



4(a + 7) J VV a + 7 

Proof. We first expand the fraction 

Z%/Zr, = where a; = (a + 7)"^/^ 

The diagonal elements of a GUE^ matrix are independent Gaussians with mean and variance 
2Af(a+7) • Using this, we get easily for any e, any N and any a, 7 > 

Zh/Z^";^ = E;v,GUE. exp{e^Mi(x)} = exp {^e' ■ (4(a + 7))"'}, 
where Eat^guEi^ denotes expectation w.r.t. P^^J^- Similarly, we get for any N and any a, 7 > 
Zn/Z^I"" = E^,GUE„ exp {7Mi(x)2} = (1 - 7/(0 + 7))"'/' • 

□ 

Definition 7. For a; > 0, the probability measure a^^ on M given by 

2 

cr^(t)dt := 2 V^^^ -*^l[-<^,t^l(a;)c^i 

is called (Wigner's) semicircle law (with parameter cj). 

By equation , Pat is a mixture of . We show first that the statement of Theorem [S] is 
true for each e G M if we replace P^ q by P^. Eventually we will use Lebesgue's dominated 
convergence theorem. 



Proposition 8. Let p^^f denote the k-th correlation function of Pf^ and set lj — y . 

a) For any e G M, any k and any continuous, bounded g : — > M we have 



lim / g dp}f = g d{a^ 



b) We have for any e and any k, 

J^oo + N^y iv^) = - '^■))i<^.^<^ 

locally uniformly for all ti, . . . t^ and uniformly for a varying in a compact subset of 
(-w,cj). 

Proof of Proposition A proof of the first part can be found in [Joh98j . For the second part 
we use orthogonal polynomials. Note that the polynomials orthogonal to a Gaussian weight 
with non-zero mean are normalized shifted Hermite polynomials. Let vr^-^'' denote the j-th 
Hermite polynomial orthonormal w.r.t. the weight e-^i'^+'y)^^ _ 

It is easy to check that the set of polynomials orthogonal w.r.t. the weight e"^*-'^"'''^^* +"^^7* 

(N)* 

i 

:= e(-"^V2A^)^(N)(^ _ ^/^/27V) with u' := ^/{a + 7) and u" := u;'^/4. (14) 



are the polynomials [T^j ) , where 
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The ensemble Pf^ is determinantal, i.e. 

...,tu) = {N- k)\/{N\) det *i))t=i ' (15) 

where Klj{t,s) = T.^~Q^f^*{t)-Kf^*{s). From (m we get 

K*^(t, s) = e('^"^')/^K7v(i - uj'e/2N, s - oj'e/2N), (16) 
where K]\f denotes the kernel corresponding to the ensemble P^^J^- Hence we have 



aUa) Na^a)' Na^a) 



g(a;"s2)/iv t-cj'ea^{a)/2 s - w'ea^(a)/2, 



= — T~r^^i"' + i^T — t^^^ + itt — tt)' (1"^) 

where t' := t — uj' £a^{a) /2 and s' := s — u' ea^{a) /2. It is well-known that 

1 / s' , sin(7r(t' - s')) , , 

hm — — i^:Ar(a+— — rT>« + T7 — rr) = — -r, tt^- (18) 

For aproof of ([T8|) see e.g. |Dei98l Chapter 8] or Theorem [26l Since limAr^oo exp {(a;"e^)/A^} = 
1, we get from (HID and ^ that 

]_ , t s ^ sin(7r(f - s')) sin(7r(t - g)) 

lim — r^-fV/v« + T7 — r^if^^ + TT — r^) = t"; 7^ — = 1 ^ — • (19) 

N^ooa^{a) iVcr^(a) Na^{a)' 7r(t' - s') 7r(t - s) ^ ' 

Now, by (fT9]) and (fT5]) . the second assertion of Proposition [8] follows. As ([T8]) is true locally 
uniformly in t' , s' and uniformly in a G /, / C [— compact, we get (jl9p locally uniformly 
in t, s and uniformly in a € /. □ 

Proof of Theorem [5l By equation (I13p and Lemma [6] we know that 

P7v(x) = / Plfix)p{e)de, (20) 
where p is an A^-independent probability measure on M. Using Fubini's Theorem, (j20p implies 

gdp%= / / g dp''^" p{e)de and p%{ti, . . . ,tk) = / p'}f{ti, . . . ,tk)pie)de, 
and hence for each compact AT C M'^ and each compact / C (— w,w) 



^1 , ifc 



sup o-aj(aj /9^(a + — — a + 



Naja)' ' A^(j,, 



_)_det (S(t.-t,))i<„<, 



(21) 



where we stick to the notation of Proposition [HI Theorem [5] will follow from Proposition [8] if 

/jgfc g dp]f and sup^g^ ^e/ Pa^^C^i, . . . , s^) , := a + U/ {Na^{a)), are uniformly bounded in 

£. The uniform boundedness of J^k 9 dp'^f is immediate as g is bounded. 

To show uniform boundedness of p'lf{si, . . . , Sk) uniformly in e, t and a, we proceed as in 
the paper by Pastur and Shcherbina |PS08j . Since all correlation functions are nonnega- 
tive, we see by Sylvester's criterion from the determinantal relations (jlSp that the matrix 
{K*j^{ti, tj))i<ij<k ='■ ^ is positive semidefinite and can hence be written as A = for some 
matrix B. Now using Hadamard's inequality we get 

k k k 

det A = (det < H '^^j'^ = 11 ^JJ- 

our case this reads 

j=i 1=1 j=i 

k k 

p'^/ih, ...,tk)<{N- ky./{N\) H KN{tj,t,) < C^Hp'j^'itj), (22) 

i=i 3=1 
where C is a constant such that C > N/ (N — k). Using p4|) . we get 

N-l 
j=0 



N-l 

^ E - ^'e/2iV)2e-^(°+^)(*-'^/2^)' = p'jf^'^-it - Je/2N), 



j=0 

1,GUE„ - 



where pj^ " is the first correlation function of the GUE^. In Proposition [19] we will prove 
that there is a constant C such that for all N and all t we have p]^^^^"it) < C. This 
estimate together with (j22p finishes the proof of Theorem [5l □ 

3. The Associated Random Matrix Ensemble 

In this section, we start with the investigation of our main model. Let h he a continuous 
even function and Q a strictly convex symmetric function and assume that 

Pk,Q{x):=-^ n |x,-x,f e-^S.-«(^^)-S-.'^(^-^^), (23) 

N,Q l<i<j<N 

defines the density of a probability measure on M^, where 

-'^^ l<i<j<N 

denotes the normalizing constant. This is for example the case if h is bounded below. 
We will frequently use the notation 

h^{s):= j h{t- s)dp{t), h^^r-= j j h{t- s)dp{t)dp{s). (24) 

for a compactly supported probability measure p on M. For the statement of the next lemma, 
M.]. will denote the set of compactly supported (Borel) probability measures on M. 
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Lemma 9. Let h : M — > M be even, twice differentiable, bounded and such that h"{t) > —aq 
for all t. Define : Ai^ — > -M.]., Th{fj.) as the equilibrium measure to the external field 
t^ Q{t) + hf,{t). 

Then has a fixed point, i.e. there exists a probability measure fig which is the equilibrium 
measure to the external field 1 1— )• Q{t) + f h{t — s)dfiQ{s). 

Proof. We will apply Schauder's Fixed Point Theorem, which states that each continuous 
mapping T : C — > C of a compact, convex and non-empty subset C of a Hausdorff topo- 
logical vector space has a fixed point. 

We consider the topological vector space A4{K) of all signed finite Borel measures on some 
compact interval K of M, equipped with the topology of vague convergence. This topology 
is metrizable and hence the space is Hausdorff (see [ST971 Chapter 0]). The subset M^{K) 
of all Borel probability measures on K is non-empty, convex and compact. The compactness 
follows from Helly's Selection Theorem. We will further restrict to measures n which are 
symmetric around 0, i.e. /x(A) = IJ,{—A) for all Borel sets A. It is easy to see that this subset 
still fulfills the assumptions of Schauder's Fixed Point Theorem. 

Now we show that since h"{t) > —aq and h is bounded, the support of the equilibrium 
measure to the external field Q{t) + h^{t) is included in a compact set which can be chosen 
to be independent of fi. Indeed, by Theorem [35} the support of the equilibrium measure for 
Q{t) + h^{t) is the smallest compact set K (w.r.t. inclusion) of positive capacity maximizing 
the functional 

K ^ FQ+h,{K) = log cap (K) - 2 j Q{t)dujK{t) -2 j h^{t)dujK{t). (25) 

= Fq{K) ~ ^ j h^{t)dujK{t), in particular we have 

FQ+/i^(supp/iQ) > FQ(supp/iQ) - 2\\h\\oo G M, since \h^\ < \\h\\oc- (26) 

As Q is convex and symmetric, supp/ig is a symmetric interval (see Theorem [35]) . Because h 
is twice differentiable, h' (and by assumption also h) are bounded on any compact set. Hence, 
if we choose a probability measure with compact support, /i^ is two times differentiable 
and {h^y = {h")^. By the condition h"{t) > —aq, Q{t) + hf^{t) is convex for each compactly 
supported fi. Theorem [351 implies that the support of the equilibrium measure to Q{t) + h^{t) 
is a symmetric interval, say \—lfi,l^- Using Lemma [301 we can rewrite (|25p for an arbitrary 
symmetric interval [—1,1] as 

FQ+h,{[-lA) = log(V2) -2 j'^Q{t)—j^==dt-2 j\{t)—j^==dt. (27) 

Since Q is strictly convex and symmetric, we have Q{t) > aqt^ + C for some C G M and (j27p 
implies (using that the variance of ^[-i^i] is l"^ /2) the inequality 

Fq+h,{[-l. I]) < log(//2) - aql' -C + 2\\h\\^, (28) 

which holds for any /j,. Comparing (j26p and (j28p . we see that 

Fq+h^isuppfiq) > Fq+h^{[-l,l]), 

for all I > L, where L > does not depend on /i. Hence such an [— /, /] can not be the support 
[— ^/i) ^/j] of the equilibrium measure for Q + /i^. Hence < L for all compactly supported ^. 
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We have thus seen that Th maps the set A4l{K) of symmetric probabihty measures supported 
in K into itself, if K is chosen large enough. It remains to show continuity of this map. Since 
we deal with a metric space, it is enough to show that by Th, converging sequences are 
mapped to converging sequences. Let {nn)n C ^A^{K) be a sequence converging vaguely, or 
equivalently, weakly to a probability measure /u. Denote Th{fJ-n) =■ ^n- Define the sequence 
of external fields Vn{t) := Q{t) + h^^{t) which converges pointwise to V{t) := Q{t) + h^{t). 
We may assume that this convergence is uniform: By Theorem 1331 the equilibrium measure 
does not depend on values of the external field outside of its support (from which we know 
a priori that it lies in a certain compact set). Since h' is bounded on this compact set by 
some constant, say C, we also have < C. This implies that the sequence of functions 

{hf_i„)n is uniformly Lipschitz and hence equicontinuous. It follows that the sequence {Vn)n 
is also equicontinuous. Since their domain is a compact and Vn converges pointwise, the 
equicontinuity implies uniform convergence by Arzela-Ascoli's Theorem. 
Since all f„ are supported on the same compact set, it follows that (fn)n is tight and hence 
has a weakly converging subsequence {i'nm)m- We will prove that this limit measure, say v' , 
is in fact v = Th{n), the measure belonging to the external field V, and does not depend on 
the particular subsequence. It follows that the sequence (fn)n converges to u weakly as weak 
convergence is metrizable. 

From the uniform convergence of Vn towards V it follows by Theorem [M] 1. that 



C7^"-(s) = / log \t - s\-Uiyn At) 



converges uniformly (on C) towards U'^{s) := Jlogjt — s| dv{t). On the other hand, by 
Theorem 1341 2. we have for all s G C except a set of zero capacity 



lim [/'^""(s) = ^/'"'(s) = \og\t- s\'^ dv'{t). 

m—>-cx) J 

Hence U'^{s) = U^' {s) almost everywhere on C. Theorem 1341 3. yields that v = v' , implying 
that the sequence (fn)n converges weakly to v. As is a continuous mapping, Schauder's 
Fixed Point Theorem yields the existence of a fixed point. □ 

Remark 10 (Uniqueness). So far we did not prove that this fixed point of is unique. 
Uniqueness will follow for the class of ensembles from Theorem [TJ For those ensembles we 
will show that the first correlation measure converges weakly to any fixed point, which shows 
uniqueness. 

We proceed by decomposing the additional interaction term. Let h be as in Lemma[9l Choose 
a fixed point fiq as in Lemma [9l We will stick to this measure from now on and write 
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h 



instead of fig. We set using the notation ([J 

h{xi - Xj) 

i<j 

N ^ 1 ^ 

j=i «j=i 

= — ^^/i/i - y/'-(0) + N 2_^h^{xj) - U{x), where 



2 



1 ^ 



Now we can rewrite q as 



pU-) = Y^ n k.-x/e-^S.-^(-^)+"(-), (30) 
^N,V,U 



^<i<j<N 

where we defined the external field 

V{t) := Q{t) + h^{t) 

and absorbed the constant exp{ — (A^^/2)/i^^ — (A^/2)/i(0)} into the new normalizing constant 
Zn,v,u- We will from now on work with this representation of the density of P^q- The 
proofs of Theorems [1] and [3] rely on comparison with the unitary invariant matrix ensemble 

PN,vix) = -^ n |x.-x,f e-^S."i^(^^). (31) 

l<i<j<A' 

We will show that in the large limit, the correlation measures in the global scaling as well 
as correlation functions in the local scaling, are the same for Ppf q and Pn,v- In this sense 
the quantity U will turn out to be negligible. 

4. Concentration of Measure Inequalities 

We will frequently use the following well-known concentration of measure inequality ( |AGZ10"t 
Section 4.4]). 

Theorem 11. Let Q be an external field on an interval I = {a,b) (possibly unbounded) with 
Q" > c > on I. Then we have for any Lipschitz function f on I and any e > 

N N 



E^,Q exp {e( ^ fixj) - En,q ^ f{xj)) } < exp {^^}, 
where for any Lipschitz function f we denote its Lipschitz constant by \ f\^ (on I). 
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N 



Remark 12. In [AGZlOj . only the case (a, 6) = M is stated. As the proof for general (a, 6) is 
completely analogous, we do not give it here. 

Theorem [11] yields a concentration inequality for linear statistics around their expectations. 
However, we rather need concentration around their "limiting expectations". It is well-known 
(see e.g. |Joh98l Theorem 2.1]) that for bounded and continuous functions 

1 ^ r 

lim ^e5^/(x,) = / fit)df,Q{t), (32) 

where fiQ denotes the equilibrium measure to Q. We need to quantify the rates of convergence 
in ([32]) . The following is a special case of a result in |Shcll| (see also [KSlOj ) . 

Proposition 13. Let Q be a convex external field on M which is real analytic in a neighbor- 
hood o/supp(^q). Let f be a function whose third derivative is bounded on a neighborhood 
of supping). Then 

N 

mm. 

where C does not depend on N or f and \\ ■ ||oo denotes the bound on the neighborhood of 
supp(/iQ). 

From Theorem [TT] and Proposition [T3] we immediately get the following concentration in- 
equality. 

Corollary 14. Let Q be a real analytic external field with Q" > c > 0. Then for any Lipschitz 
function f whose third derivative is bounded on a neighborhood o/supp(/iQ), we have for any 
e > 

EN,Qexp{e{J2f{xj)-N J f{t)dfiQ{t))}<exp{'-^ + eC{\\f\\oo + \\f'^'^\\oo)}. 

Remark 15. Proposition 1131 and Corollarv 1141 remain true up to an error of order e"*^^ if we 
replace M by an interval / which covers the domain of the equilibrium measure /ig . It is well- 
known (see e.g. jPS08] . p3Gllj ) that changing the external field outside a small neighborhood 
of the equilibrium measure results in a change of the first correlation function of order e~'^^ 
for some c > 0. We will prove this in Lemma [28] provided that / is large enough. 

The next lemma gives, using Fourier techniques, a representation of the bivariate statistic l/l 
in terms of certain linear statistics. A similar idea is used in |LP08| . 

Lemma 16. The following holds 

lA{x) = -j= I UN(t,x) h{t)dt, where 



2V27r 

N ^ N 



UN{t,x) :='S^cos{tXj) — N f cos{ts)diJ,{s) + y/^'S~]sm{tXj), h{t) := / e ^^'^h{s)ds. 

frf J TZf v27r Jm. 
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Proof. Recall from ([29]) that 



1 ^ 

^{x) = ~9 ( X] ^(^« ~ ^i) ~ + ~ Vm] )• Note that 



Xj Xfc) - 

Atx 



2v/27r 



J{Xj-Xk)t 



h{t)dt = ^^=J \uNit,x)\'^ h{t)dt, 



with UN{t,x) := e**^J . Writing U]\f{t,x) := UN{t,x) — N f e**''(i/u(s), it is not hard to 



check that 



U{x) 



1 



2\/27r 



UNit,x) h{t)dt. 



(33) 
□ 



Note that we can write 

With the help of the representation ()33p we shall bound this ratio of normalizing constants. 

Proposition 17. If the constant aq is large enough, then there exist constants Ci, C2 > 
such that for all N 

0<Ci< Zn,v,u/Zn,v = IE7V,vexp {U{x)} < C2. 
Proof. We start with proving the lower bound. By Jensen's inequality we see 

E7v,v exp {Z//(x)} > exp {ENyU{x)}. 

Using Lemma[16]we show that the expectation of U is bounded in A^. Fubini's Theorem gives 
1 




Ny 



Uj\f{t, x] 



h{t)dt 



cos(txj) — / cos{ts)dfi{s)\ + Eat vl^^ sin(txj)p 1 h{t)dt. 



By Corollary [T^ the terms in the parenthesis are bounded by a polynomial function in t, as 
|cos(t-)|£ , |sin(t-)|£ < t and ||cos(t-)^^^ ||oo, ||sin(t-)^^^ ||oo < Ct^. Hence, h being a Schwartz 
function, we have '&i\fylA{x) > —C for some C > 0. Thus the lower bound follows choosing 
C7i :=exp(-C7'). 

For the upper bound we will again use the representation of Lemma [T6l Recall that since 
h is even, h is real-valued. Define h+{y) := max{0, and h-{y) := max{0, — such 
that h = — h-. For h- = 0, which corresponds to the case of a positive definite h, there 
is nothing to prove, so assume that /i_ / 0. 
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^ 1 /2 

Introducing H- := {h ) > 0, we obtain by Jensen's inequality and Tonelli's Theorem 
EAr,yexp| - (2\/2^)"^y" h{t)\uNit,x)\'^dt^ <EN,vexp^{2V2^)-^ j H^{tf \uN{t,x)\^ dt^ 

< y (/7_(t)/||F_||^i)Ejvyexp{(2\/2^)-^||i7_||^ii/_(t)|njv(t,x)|^}(it. (34) 

Abbreviating '■= (2-v/27r)^^ and using the Cauchy-Schwarz inequahty and repre- 
sentation (j33j) . we find 

En,v exp {KhH^{t)\uN{t,x)\^] (35) 

N r 2 

< E]ly exp ^2KhH^{t)^ ^ cos(tej) - iV / cos{ts)dnis) | (36) 

^ 2 

X E^y/'^, exp {2KhH^{t)\Y^ sin(te, ) | } . (37) 



Since by Corollary [HI the distributions of Ylf=i cos{tXj) — N f cos{ts)dfi{s) and J2f=i sm{txj) 
are sub-Gaussian, we obtain for example for the first term for any e > 0, 



N 



ENyexp^e■^/2J{iJ^^t){^cos{tXj)-N / cos{ts)dfi{s))^ 

< exp {e^ . 2KhH^{t)t\2av)'^ + eyj2KhH^{t)C{l + t=^)}, (38) 

where ay := min^ V"{t) > 0, C does not depend on t or A^. For aq large enough (hence 
av large enough), we have 2KhH.{t)t^{2av)~^ < 1/4 for ah t. Since F_(t) = /^^/^(t) is 
decaying rapidly, ^J2KhH^ {t)C{l A-t^) is bounded in t. Summarizing, if uq is large enough, 
we can bound (p8|) by 

exp{ce^ + eC} 

with < c < 1/4 and c, C do not depend on N or t. We conclude that (p6|) and (f37|) and 
hence ([35]) are bounded in N . Finally, since /i is a Schwartz function, it follows from ()34p 
that 

E7v,\/exp{— y /i(t)|ti7v(i, 2;)p(it} < C 

for some constant C > independent of N . This proves the upper bound and hence the 
proposition. □ 

Remark 18. The proof of Proposition [T7] actually shows that for each A > there is a 
threshold a^{\) > and constants Ci,C2 (depending on A and a^) such that 

< Ci < E^y exp{A^/(x)} < Ca, if ag > a^(A). 



15 

5. Bounding the First Correlation Function 

The first aim of this section is to prove that p\; qit) < C uniformly in t and N. We wih 
eventually prove the following proposition which can be found in |PS08llPS97j . We use in the 
proof arguments and results from [Joh98[ IPSOSt IPS971 IDei98j . For the sake of completeness 
we shall include them here. 

Proposition 19. Assume that Q is three times continuously differentiable and satisfies 
Q{t) > (1 + (5)log(l + i^) for some 6 > and all t large enough. Then there is an N- 
independent constant C such that for all t and all N 

Before stating the next concentration inequality, we need to introduce some notations. 

Kn,q{x):= Y1 kQ{xi,Xj), kQ{t,s):=\og\t-s\-^ + ]^Q{t) + ]^Q{s), (39) 

l<i^j<Ar 

Fq:=Iq{p), VqW := -log(t^ + 1), where /q(/x) is defined in ([92]). (40) 

From the simple inequality |t — s| < + iV + 1 we conclude log \ t — s\^^ > log(t^ + 
l)(s^ + 1) and hence 

kgit, s) > (l/2)V'Q(t) + (1/2)V'q(s). (41) 
We also note that since V is an external field, there is a constant eg such that 

Mi) > CQ. (42) 
We define a generalized unitary invariant ensemble on (or some compact [a, b]^) via 

Pi^g,,(x) := n \x.-x/e-'^^U<^^^^^^Uf(^^\ (43) 

'^N,QJ l<i<j<N 



where iV, M E N and / is a continuous function with \f{t)\ < Q{t) for t large enough. Usually 
we have M = iV or M = iV - 1. If M = A^, we will write Pn,qj instead of P^\qj- If / = 0, 
we write Pffq- The following result is due to Johansson. 

Proposition 20. Let 

1 



AN,e ■■= ja; E : -^Kn,q{x) < Fq + e 

Then there is some constant C such that, if\miN^ooN/MN — >• 1, 

P^\q{^^ \ An^s+o) < Ce-"^^ for all N > No{e) and all a > 0. 

Proof See |Joh98l Lemma 4.2]. □ 

We abbreviate in the following pj^q = /ojy • The first step in proving Proposition [19] is to 
determine the decay of p\;{t) for t large. Fortunately the precise decay is well-known in the 
theory of invariant ensembles. The following lemma can be found in several works including 
[Joh98[|PS08llDei98j . We follow |Joh98| . 
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Lemma 21. Let Q be a continuous function satisfying Q{t) > (1 + (^)log(l + t^) for some 
5 > and all t large enough. Then there is a constant C > such that for all t 

Proof. We will from now on drop the subscript Q, defining 

(x) := — ^ Y\ ~ Xjl^ e~*^^J=i '^'•^^^ and abbreviating pjy •= Pn,q^ compute 

l<i<j<N 



j=i J 

(44) 



Since adding a constant to Q does not change the ensemble, we will assume that Q > 0, 
which corresponds to considering the potential Q + Cq, where Cq denotes a lower bound of 
Q. Setting Z := J e~^^^^dt we get by Jensen's inequality 



7V-1 , ^ N-1 



I exp{2 \og\xj-t\ -NQ{t)}dt > Z exp / log|xj -t| - {N - l)Q{t))e-^^^Ut} . 

i=i j=i 

Since Q > 0, we get 

\og\t-Xj\e~'^'-*Ut> / log\t-Xj\e~'^^^Ut> / log\t - Xj\ dt =: Ci. 

J Xj — l J Xj — 1 

Summarizing we have 

Z^/Z^_^ > Zexp{-2(iV - l)Ci}. (45) 
Using the inequality {xj — t)"^ < (1 + a^|)(l + t^) gives 

iV-l N~l 

n i-^-'f) ^ (i+tTK-'i n (46) 

j=l j=l 

As before, we can assume (otherwise we add a constant) that Q satisfies Q(t) > (l + S) log(l + 
t^) for all t and some 6 > 0. Using notations (j39ll40p and inequality (jSJ, this condition yields 

N-l 

Kn-i,q{x) > 6{N - 1) 5^ log(l + xj f. 

i=i 

Proposition [20] shows that for A large enough we have 

N-l 

Pn-i,q{ E + ^i)' ^ ^ ^'iO'-i,Q(i^iV-i,Q(x) > M(iV - 1)N) < e-'^^^' (47) 
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for some constant c > 0. Prom this we conclude that for A large enough 

N-l N-1 

K~'{ n (1 + ^^')) ^ ^^"^ + + ^i)ln"- (i+.^)>e^-)- (SZI) gives that 

+ ^^■)' - '^^(^ - ^ ^ exp{-iV2 _ ^y^y 

From this bound it is easy to see that E^-i ( Y\^-,\l + x'^Mj,. , ,N , ) is of order 
exp{— A^^}. Hence we have 

^^n'^{ n + ^i)) - exp{c^A^} for some c. (48) 

In view of (03]) we find combining (gS]), ([3S]) and (|3S]) 

pjv,Q(i) < exp{CiV} exp{-iV [Q(t) - log(l + t%}. 

□ 

From the previous lemma we easily deduce the following important corollary (cf. |Joh981 
iPSOSllDeiM] ). 

Corollary 22. Let Q be as in Lemma [2ll Then there are L,C > such that for all t with 
t > L, we have 

p]^it) < exp{-CNQ{t)}. 
Furthermore we need the following lemma from | Joh98| : 

Lemma 23. Let Q be as in Proposition\19l Let t S [—L, L], where L was defined in Corollary 
\22l Let < (5 < A^^^. Then there exists a constant C > such that for all \s\ < 6 

pUt)<CpUt + s) + 0{e-''). 

Proof Writing := [-L,L]^"^ we split up (jM]) as 

N~l 

pUt) ={Z^^i/Z^)E^-'{ n {xj - t)'lA,(x)) exp{-ArQ(t)} 

i=i 

N-l 

+{Z^_,/Z^)E^-\ n {xj-t)Wc^{x))eM-NQ{t)}. (49) 

Our first aim is to show that {Z^_i/ Zj^)E^-^ {Y\fj^^ {xj - tflAiix)) exp{-NQ{t)} can 
be neglegted as A^ — )■ 00. To see this, we repeat the considerations leading to (|18]) and 



find combining with (I45p that for some A large enough we have for some positive constants 

N-l 

^Ej^-i( n - tf (x))e-^«W < e(^i+^)^P^-nAi)e^'°s(i+^') + e'^^^^. (50) 
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Here we also used that \t\ < L and Q > 0. Note that Ci does not depend on L. Since x S A^, 
we have for some j, \xj\ > L and hence by CoroUary [22] we can find L large enough such that 
for some C4, c > 



'{Al)<N [ p]^{t)dt<exp{-Nlog{l + L^)-{Ci+A)N-C^N}. (51) 

J\t\>L 

Inserting (fSTj) into (f50|) yields 

i=i 

for some positive constant C. Hence we can write 

N-1 

pUt) <^^K-\ n -t)'lA.(-))e-^«(ne-^^ (52) 
= fE_L,-NQit) f Yl {x, - tf n (^^ - x,)^e-S.- + e"^^. (53) 



From the Lipschitz continuity of Q on [— -L, L] follows that 

sup \NQ{t) - NQ{s)\<C5N. (54) 

|t|<2-L,|s|<<5 

Now we can make the change of variables Xj = yj — s in ([53|) and use ([5^ to get 

p]v(i) < exp{C6N^}plfit + s) + exp{-CiV}. 

Since by assumption < 6 < N~'^, we get the desired estimate. □ 

Before bounding p]y also for \t\ < L, we need one more lemma due to |PS97j . 

Lemma 24. Let he the reproducing kernel for the external field Q, i.e. K]\f(t,s) := 
'Yli=o' where ip^^^ are the orthonormal functions w.r.t. the weight e~^^ . 
Then there are constants C, C" > such that 

j {t- sfKN{t,sfdtds <C 

and for a = 1,2 

1 1 it-srKMit,sfds\ < C7'«^Zi(i)2 + <^(f)2). 
Proof. It can be checked that using the well-known Christoffel-Darboux formula, we get 

J\t- sfKNit,sfdtds = 2(j^_i)', 

where jf^ := J tipi^\t)il)''^\{t)dt. Furthermore we have for \t\ > L hy Corollary [22] 

^Pritf<Nplj{t)<Ne-^'"^^'\ 



19 

which imphes | Jjv_i| ^ C for some C. This proves the first part of Lemma Y2M. The second 
part can be proved similarly using 



□ 



Now we prove the uniform bound for pjy (also due to |PS97] ). 
Proof of Proposition The starting point is the simple identity 

dp]^{t) _ dp]^{t + s) 



(55) 

dt ds ^ ' 

s=0 

Differentiability of can be seen easily using its determinantal representation. Changing 
variables to Xj — s for all j, we can write 

p],{t + s) = -^ [ e-^Q(*+^) n (x, - x,f n e-^0(^+-^) {t - x,f dx. 

Differentiating w.r.t. s and setting s = later, we get 

"^^^^^^ - NQ\t)pj,{t) - NiN - 1) / Q'is)p%it, s)ds 



dt 

= -Q'{t)KN{t,t)- J Q'{s) {KN{t,t)KN{s,s) - KN{t,sf) ds. (56) 

Here we used the determinantal representation of p^. Integrating this relation w.r.t t we get 
using Tonelli's Theorem 



Q'{t)KN{t,t)dt = 

and hence ([56]) reduces to 

dp],it) 



dt 



{Q'{s)-Q'{t))KN{s,tfds, (57) 



where we used Kiy{t,t) = f Kiy{t, s)'^ds. We split integral (f57|l into two parts corresponding 
to |s| > L and |s| < L, where L was defined in Corollarv I22i Now we observe that from 
< N{N -l)p%{t,s) = KNit,t)KN{s,s) - KN{t,sf the simple inequality 

KN{t,sf < KN{t,t)KNis,s) 

follows. Hence we can estimate 



\[ {Q\s)-Q\t))KNis,tfds\<[ \Q'{s)\KNis,tfds + \Q'{t)\[ KN{s,tfds 

J\s\>L J\s\>L J\s\>L 

<KN{t,t) \Q'{s)\KN{s,s)ds + \Q'{t)\KN{t,t) KN{s,s)ds. 

J\s\>L J\s\>L 

Now from Lemma [23l we can conclude that there are constants C, C" > such that K]\f{t, t) < 
CN'^' for all t. This can be seen easily as follows. For |t| > L we can use Lemma [2T] and for 
\t\ < L, Lemma [23] in combination with J KN{t,t)dt = N and -R'at > gives that there are 



20 



constants C, C" > such that KN{t,t) < CN^' . Using again Lemma [2T] we conclude that 
there is a positive constant C such that 

\! {Q'{s)-Q'{t))KN{s,tfds\<C. (58) 

J\s\>L 



To estimate the second part of the integral (j57p . we use Taylor's expansion centered at t 
yielding 

Q'{s) - Q'{t) = {s- t)Q"{t) + ^{s- tfQ"'{0, 
for some ^ between s and t. Bounding Q'" and Q" on [—L, L] by C, we get 

\[ {Q'is)-Q'{t)) KNis,t)'^ds\<\Q"{t)\ [ \s-t\KNis,tfds + ^ [ \s - tl"^ KN{s,tfds 

J\s\<L J\s\<L 2 J\s\<L 

< C'«!,(t)2 + V^j,-)(t)2) +C7"«l,(t)2 + <'(t)2), (59) 

where we used Lemma [Ml Combining (j59p with (j58p we find using (j57p that there are positive 
constants C, C such that 

dpUt)/dt < C{i^^^l,{tf + <)(t)2) + (60) 
To infer from this bound the boundedness of p]^, we use for a < b the simple equality 

Plrii) = [ PNit) - Phis)ds + / pl,{s)ds. (61) 



Since p% > and / p]^{s)ds = 1, we conclude that the second integral in ([611) is bounded by 

(a — b)^^. For the first interval we can write p\f{t) — p\f{s) = -^-^;^du — /J -^^i^du. We 
can assume without loss of generality that \t\ , |s| < L (otherwise we get a bound by Lemma 
[2T]) and hence using (j60p in combination with the orthonormality of '0''^' gives the desired 
bound. □ 

We finish the section with a useful bound on the first correlation function pjy q / °f unitary 
invariant ensemble Pn,qj (see (03])). 

Corollary 25. Let f be continuous and bounded. Then there is a constant C > 0, only 
depending on Q such that 

Proof. We use the identity 

^ e-NQ+f 
Pn,qA^) = NXN{e-^Q+f,ty 

where X]\[{e~^'^^^ , ■) is the so-called A^-th Christoffel function to the weight e~^^~^^ (see 
[TotOO] for references and more information on Christoffel functions) 

\N{W,t):= inf [ \PN^i{s)\^W{s)ds, (63) 

PN-l{t) = lJ 

where the infimum is taken over all polynomials Pn-i of at most degree — 1 with the 
property that P/v-i(t) = 1 and W denotes a weight function on M. It is obvious from ([63l) 
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that XnIWi, •) < Xn(W2, •) if Wi < W2. Then the corollary follows easily by e'^^'^^-^^^^ < 
g-NQ+f < g-iVQ+||/||oo and Proposition [m □ 

6. Proofs of Theorems [U and [3] 

We first cite a general result by Levin and Lubinsky ([LLOSI, Theorem 1.1]) about bulk 
universality for unitary invariant ensembles. Recall the definition of p% q j following (j43p . 

Theorem 26. Let Q he a continuous external field on the set S C M, which is assumed to 
consist of at most finitely many intervals. Let f he a hounded continuous function on S. Let 
Kn denote the kernel 

N-l 

i^^(t,.) = 5^^f^(t)^f (.), 

j=Q 



where y^j^j o'^e the orthonormal functions to the weight e . Let J he a closed 

interval lying inside the support of fj,Q. Assume that fig is ahsolutely continuous in a neigh- 
horhood of J and that Q' and the density fig are continuous in that neighhorhood, while 
fiQ > there. Then uniformly for a G J and t, s in compacts of the real line, we have 



^.^ + KN(a,a) ' " + KNla,a) ) _ sin (7r(t - g)) 

iV->oo Kjq (a, a) Tr{t — s) 



We use a notion of bulk universality which slightly differs from (j64p , namely we scale by the 
limiting density instead of using the A^-particle density. The following obvious corollary 
is a translation of Theorem 1261 into this setup. 

Corollary 27. Let Q, f and hq he as in Theorem\2(A Then hulk universality as defined in 
([2]) holds for the unitary invariant ensemhle Pn.qj- 

Proof. The corollary follows from the well-known determinantal relations for unitary invariant 
ensembles, the local uniformness of the limit (j64p in t, s and the fact that by [TotOOl Theorem 
1.2] we have uniformly in compact proper subsets of supp/iQ 



lim ^KN{a,a) = lim PN,Qj{a) = pqia). 



□ 

We will prove Theorems [T] and [3] together by comparing the correlation functions of the 
ensembles q (see ([30]) ) and Pn,v (see (j3T|) ). We start with p^y, the fc-th correlation 
function of Pn,v- We obtain p%y{a + jy^^, . . . , a + jy^^) ^ /c-marginal, integrating the 
density 

PN,v(a + T7-TT' • • • + . ,Xk+l, . . .,xn) 

Nfi{a) Np{a) ' 

over Xfc+i, . . . , xat. We have k fixed eigenvalues at positions a + ]y^^, . . . , o + ^y^^ and 



N — k random eigenvalues. We first rewrite p%y ( a + jf^j^, • • • , a + j^^i^ j in terms of these 
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N — k random eigenvalues as follows. 

N 



exp{-^En«+^)+2 E log|^|} (65) 
exp{2 ^ log |a + ^^'^^^^ - Xj\]dxk+i ■ ■ ■ cLxn 



i<k, j>k 



(r, 1 r:.A — TT It- - r-|2p-^^f=fe+i^(^i) 



= F{a,t) ;^~^'^ Ej^_fc yexp{2 V log |a + — ^ - x^l }, where (66) 

F(M):=exp{-A'tl'{» + ^)+2 l„g|^|} (07) 

is the factor ([65]) . which depends only on the fixed particles, and 

1 -TT , ,9 -ArV^JV 

Af-fc,V k^\<i<j<N 

As before, the subscript N — k indicates that P^_k y is a probability measure in — A; 
variables, whereas the superscript indicates that the factor in front of the external field 
term Vixj) of P^_k y is A and not N — k. We keep the labelling Xk+i, . . . , xjy. 

Setting 

L^_,y{a,t,x):=2 J] bg |a + — ^ - | + log [F(a, t)^^] , (68) 
^<t:1>k ^^^''> ^^'^ 

we get from ([66]) the equality 

Pn,v (q + jy^^^^^ , ■ ■ ■ , a + ^^^^^^ ) = Ej^_;. yexp{L]^_fcy(a,t,x)}. (69) 

Similar as in (j66|) , we see that the fc-th correlation function P^'^q of q at a + j^^^ , . . . , a + 
s can be written as 



E^_j,y exp {U{t, x) + L%_f,y{a, t,x)}, (70) 



where we abbreviated U{a + jy^) ■■■ ,a + -j^^,Xk+i, ■ ■ ■ ,xn) by U{t,x). 

In the following we shall abbreviate (ii, . . . , tk, Xk+i, ■ ■ ■ ,xiy) by {t, x) and by {t, x)j we will 

denote the j'-th component of the vector (t, x). Furthermore, for the sake of brevity, we set 

Ra:= L^_k,v{a,t,x) and R := L%^,,y{0,Np{0)t,x). (71) 
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Note that R arises in the global scaling, whereas Ra appears in the local scaling. It will later 
turn out to be convenient that all the Xj's lie in a compact set. To this end we formulate the 
following truncation lemma. This procedure is well-known for invariant ensembles (see for 
instance |,Toh98j or |BdMPS95] ). 

Lemma 28. For ag large enough, the following holds: For each k there are L,C > such 
that for all N and for all ti, . . . ,tk 



p^^^iti, . . . ,tk) - f . E^_ky,L(iw{^it,x) + Rl} 



where Kfjy^ denotes expectation w.r.t. the ensemble Pffyi^ obtained by normalizing the 
ensemble P^y restricted to [— L,L]^ and Rl is the analog of R in which all integrations 
overM have been replaced by integrations over [—L,L]. Furthermore, for any external field Q 
on R, the following holds: For each k there are L,C > such that for all N and all ti, . . . ,tk 



PN,Q{tli ■ ■ ■ itk) — Pn,Q,L'{'^1^ ■ ■ ■ ,tk) 

where q ^, is the k-th correlation function of the ensemble Pn,q,l' obtained by normalizing 
the ensemble Pn,q restricted to [—L',L 



nN 



Proof. We will use the representation (|70p and show that the restriction of integrals to 
[— L,L]^ C respectively [—L,L]'^~^ C M^~'^ results in an asymptotically negligible 
error. For E,j\[ye^ we use Holder's inequality to estimate 

En,v (exp{Z^(x)}l([_i,i]iV)c(x)) < (EAr,yexp{(l + (PiV,y (([-L, L]^)') ) 

where 1/(1 + e) -\- 1/e' = 1 and e > is fixed. Now E^ye^^^^^^^^^ is uniformly bounded in 
N by Proposition [T7] provided that ag is large enough. Furthermore, by Corollary [22] we get 
for the L defined there 

Pr,,v ( {[-L,Lry) <N [ p%y{t)dt <N [ e'^^^Wdt < e"^'^ (72) 

^ ^ J\t\>L ' J\t\>L 

for some C > 0. In fact, C can be chosen arbitrarily large by increasing L. We conclude 
that 



Ejv,y (exp{^/(x)}l([_i^i]iv)c(x)j < exp{-C7"iV}, 

for some C" > 0, if L is large enough. It follows by (j72p as well that the exchange of the 
normalizing constants Z^^iy and Z^_f^ y by their counterparts Zn,v,l and Z^_j^ y ^ and hence 
also the exchange of R by Rl is asymptotically negligible. 

In order to bound E]^_^ y ^exp{Z^(t, x) + R^1(^^_i^ i^^n-^'={x)^, first use Holder's inequality as 
above. It remains to estimate E]^_^ y exp {(1 + e)U{t, x) + (1 + e)i?} for some fixed e > 0. 
Again by Holder's inequality we reduce this to bounding E]^_^ ^ exp | (1 + e')Z^(t, x)| and 
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E 



y exp {(1 + e")R} for some e',e" > 0. Recall from ([33]) that 



U{x) 



1 



2\/27r 

N 



un{s,x) 



h{s)ds, where 



N ^ N 

un{s,x) = ^^cos(sxj) — N cos{st)diJ,{t) + \/— 1 sm(sx 



For any a and any ti, . . . ,tk we get 



Kit,x) < 



1 



+- 



2V27r 
1 



<- 



2V27r 
1 



27r 
1 



= fc + l J = l 

sin(sxj) + sm(stj)| h^{s)ds 
=k+l j=l 

cos{sxj) — {N — k) / cos(sii)d^(ii)| h^{s)ds 



j=k+i 

N 



sm(5Xj)| /i_(s)ds+— ^= y h-{s)ds^ 



(73) 



i=A:+l 



where we used the inequalities + < 2(a^ + 6^) and |cos| , |sin| < 1. From this we conclude 
as in the proof of Proposition [T71 that Ej^_^yexp{(l + e')L({t,x)^ < C provided that ag 
is large enough (which does not depend on k) and C does not depend on ti,. . . ,tk or A^. 
To see that Theorem [IT] also applies for Pffuy, is obvious, for Proposition \T3\ we use that 



N 



P. 



N-k 



N-k,V ~ ^ N-k,VJ 



with f{t) := kV{t) and the notation introduced in (P3|) . Proposition [T3l 
is proved in [Shell] also for the case of Pn,qj for real-analytic Q and /, hence it can be 
applied as in the proof of Proposition [T71 We may now bound IEj^_^ y exp {(1 + £")R} as in 
the arguments following ([I6]) . Recall that 



R := 2 ^ log |ti — Xjl + log 

i<k, j>k 



F{0,Nfi{0)t 



^N~k,V 

Zn,v 



where F{a, t) was defined in (j67p . Using the same Jensen type trick as in the proof of Lemma 
El we find that y/Z^y < exp{CkN} for some C. As in (|46p we get 

E^_^,yexp{{2 + 2e") J] log|t,-x,|} 

i<k, j>k 

< exp { (iV - A;) (1 + e") log (l + t?) jEj^.^^v- exp { (1 + e") ^ log(l + x|) } . (74) 



i<k 



j>k 
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Analogously to ([IHD we conclude that Ej^_^ y exp {(2 + 2e") Ej>fc log(l + x^)} < exp{ciV} 
for some c > 0. Using (|42p . it is straightforward to bound 



exp {(iV -k){2 + 2e") J^log (l + tf) + log [F{Q,N ^i{mz^-ky / Zn,v] ] 

i<k 

k 

< exp{-ciA^ - C2log(l + t^)] + CkN}, (75) 

i=l 

where ci,C2 are absolute positive constants. Since V is strictly convex, this yields 
E^_i.y exp { (1 + e")R} < e^^ and hence 
^%-ky exp { (1 + e)U{t, x) + (1 + e)R] < e^'^ 

for some C, C". From ()72p . we get that for L and N large enough 

Ej^_,,^ (exp{Z^(t,x) + i?}l([_^,^]iv)c(x)) < e-^"^ 

for some C" > and all ti, . . . , t^. 

From ([70|) . (j73|) and ([7¥|) we also obtain similar as in Lemma [2T] 

P%{ti, ...,tk)< exp{CN - ciNj2[yiti) - C2 log(l + t?)]} 

■t=i 

for some positive C, ci, C2. As before, this implies that we can assume all ti, . . . , to lie in 
some compact set. 

The second assertion of the lemma follows analogously from (j72p . (j75p and ()74p with e" = 
0. □ 

Proof of Theorems CI anc? 0. We first outline the main idea of the proof. Recall from ()29p 
that 

AT 

U{x) = -(1/2) ( ^ h{xi - Xj) - [h^{xi) + h^{xj) - h^^]). 

Assume for a moment that —h/2 is positive semi-definite, or in other words, the covariance 
function of a centered stationary Gaussian process (G'j)jg[_^^^], i.e. —h{t — s)/2 = E(GjGs). 

We may linearize the bivariate statistic —(1/2) h{xi — xj) via 

N N 

exp{-(l/2) h{x,-Xj)} = Eexp{J2G^^}, 

i,j=l j=l 

where E denotes expectation w.r.t. the underlying probability measure. By definition we 
conclude that 

N 

exp{U{x)} = Eexp {^G^^ - N / G.d^}, (76) 

provided that G. is a.s. integrable w.r.t. ^. Since we would like to apply Corollary [U] to 
the linear statistic in (j76p . we need that G. is sufficiently smooth with probability one. To 
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see this, we use the weh-known Karhunen-Loeve expansion of G. By a classical result due to 
Mercer the covariance function h admits an expansion, converging uniformly on [— L,L], 

oo 

-h{t - s)/2 = >^Mt)'W), (77) 

i=l 

where {9i)i denotes an orthonormal system of eigenfunctions of the integral kernel h with real 
and positive eigenvalues (Xi)i, i.e. 

J -{l/2)h{t - s)ei{s)ds = XiOiit) Vi. 

The Karhunen-Loeve expansion of G is then given by 

oo 

Gt = Y.X]'''^A{t). (78) 

i=l 

where (^j)j, '■= {Xi)~^^'^ j^j^9i{t)Gtdt, are independent standard normal variables. The 
convergence in (I78p is a.s. uniform on the compact interval [— L,L], see [AT07t Theorem 
3.1.2]. The a.s. continuity of Gt used for this theorem is easily checked using any of the 
well-known conditions. Since h is analytic on some domain containing the compact set, 
say A := [— L, L] x [—5, 5] C C, 5 > 0, its eigenfunctions (with nonzero eigenvalues) are 
analytic on A. Hence the uniform convergence in (j78p implies that G^, if G vl is analytic 
with probability one. Furthermore, recall that the derivative process {G'^)t^^_i^i-\ of G is 
a centered (real-valued) Gaussian process with covariance function h" /2. See e.g. jAdl8H 
Theorem 2.2.2]. 

To summarize, if —h is positive semi-definite, U admits the linearization (j76p in terms of 
linear statistics with random test functions which fulfill the prerequisites of Corollary [HI 
However, as —h is in general not positive semi-definite, we may extend the previous case by 
means of the following argument. Recall the decomposition of h into nonnegative functions 

h = (/i)+ — {h)-. By setting h'^ := (/i)+, h~ := (h)-, we get a decomposition h = h'^ — of 
h into positive semi-definite, real-analytic functions. Define for a complex parameter z £ C 

N 

1 ^ 

+ 2 ( E ^^(^^ - - [^/^(^*) + ^/^(^i) - ^mm] )• (80) 

Note that U-i = lA. We shall use the following result, known as Vitali's Convergence Theo- 
rem, which can be found e.g. in [Tit39j . 

Theorem 29 (Vitali's Convergence Theorem). Let fn{z) he a sequence of analytic functions 
on a region D C C with \ fn{z)\ < M for all n and all z £ D. Assume that liuin^oo fniz) 
exists for a set of z having a limit point in D. Then \\m.n^ao fn{z) exists for all z in the 
interior of D and the limit is an analytic function in z. 
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Let fc E N be fixed. We want to show that 

h,k / , ^1 



(81) 



tends to zero uniformly in a and t in compact sets (a inside the interior of the support of /x), 
where 



sin {nitj -tj)) 

TT{ti - tj) 



l<i,j<k 



We will first show pointwise convergence of ()8ip to and later deduce the uniformness from 
Arzela-Ascoli's Theorem. By the boundedness of Etv,^^^ (Proposition I17p and Lemma [ 
()8ip converges to zero if and only if 



EU h,k I , i 
N,v,Le PN^Q[a + 



, a + 



tk 



EN,v,Le^SHt) 



Nix{a)' 



tends to 0, where the L > comes from Lemma [28l But this means, using (I69p . (I70p . the 
abbreviation i?a,Li which denotes a version of Ra which is truncated to [— L, L] (see (j7T|) ) and 
Lemma [281 that for z = —1 



N 

N-k,V, 



■L exp {U,{t, x) + Ra,L] - ^Ny,L exp {U,]^^{t) ^ 



as — >• oo. We will prove this by applying Vitali's Convergence Theorem to the sequence 
(in N) of the following analytic functions of z: 



WN,z{a^t) ■■=K-ky,L^^v{Uz{t,x) + Ra,L] -^Ny,Le^v{Uz]^\t). 



(82) 



Introduce the domain D := {z = x + iy ^ C : x,y £ x < C{aQ)}, where C{aQ) > is a 
sufficiently small constant such that the following quantity is bounded by some constant C: 

ENyLexjp{Uc(aQ)} < C 

(the existence of such constants follows from the proof of Proposition \T7\i . First we shall 
show uniform boundedness of Wi\f^z{0',t) for all N,a,t and z £ D. By the definition of 
Uz in ([79]) and the positivity of ([80]) and ([79]) for positive z (being variances of Gaussian 
random variables) it is clear that it suffices to bound W]\f^ziO','t) for real, positive z, since for 
negative real parts of z the boundedness of W^Ar,z(a, t) is obvious. Hence we restrict ourselves 
to < z < C{aQ) only. Let and G~ denote two independent, centered and stationary 
Gaussian processes on a probability space {Q,A,P) indexed by A := [—L,L] x [— e,e] C C 
with covariance functions {z/2)h'^ and /i~/2, respectively, where and h~ are analytic on 
A. Writing Gt = Gf — J G^dfi + G^ — f GTdji and denoting by E the expectation w.r.t. P, 
we can rewrite 



E 



.fcy^i exp {Uz{t, x) + Ra,L] - ^N,v,L exp {Uz]^^{t) 



N 



=E 



E 



■ky,L { G(^t,x), + Ra,L } 



N 

ENyLexp{Y^G,^}S\t) 



(83) 
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By ([66l) we find that similarly as in (|70|) . 

{¥.N,v,L exp G^^ })" exp + ^'^^l] = pNy,G.,L{^ + ^ 



iV/x(a) 



, a + 



iV;u(a) 



(84) 

where Pn,v,g.,l denotes the determinantal ensemble on [—L, L]^ with external field exp | — 
NV{t) + Gt + Gt}. 

Fix compact sets E CM.'^ and / C supp/i°. We have 

AT Af 



sup 

teE,aei 



E 



(85) 



K-k,v,L exp { + Ra,L} - ^N,v,L exp {Y, G., }S'^(i) 

Af Af 

K-k,v,L exp + Ra,L} - ^N,v,L exp {Y,G., }S''{t) . 

j=i i=i 

Since ([8^ converges by Theorem [26] to S'^(t) locally uniformly and the term K^y^L exp {XljLi } 
is bounded in by Corollary [T3] and bounded away from by Proposition [13] and Lemma 
we see that the term 



<E sup 

teE,aeI 



N 



N 



K-kyi exp {Y,Git,x), + Ra,L} - Kn,v,l exp {^^ G,, p^t) 



(86) 



sup 

teE,aei 

converges to a.s. w.r.t. P. To show convergence of (|85p to 0, it remains to show that (j86p 
is uniformly integrable w.r.t. P. We first consider the term E^ry^iexp {X^jLi ^xj}- In view 
of Corollary UM. we need to determine the distribution of the Lipschitz constant of G^ + G^ 
and of 

||G+ + G-||oo + ||(G+ + G-)(3)||^ (87) 

on [—L, L]. The derivative processes (G"*")' and {G~y are Gaussian with covariance functions 

— {z/2){h'^)" and —{h~)"/2, respectively. Furthermore, it is well-known that supjg[_^^^]|G^| 
and supjg[_^^^]|G^| are sub-Gaussian with certain means and variances —{z/2){h~^)"{0) and 

— {h^)"{0)/2, respectively. By the same argument, ||G+ + G^||oo and ||(G+ + G^)('^)||oo are 
sub-Gaussian with certain means and the variances given in terms of derivatives of {h'^) and 
{h~). For a reference, see e.g. |AT07[ Theorem 2.1.1]. From the sub-Gaussianity of these 
quantities and Corollary [HI it is easy to see that 

Af 

EAf,y,Lexp{^G^J, (88) 

has a P-integrable dominating function, provided that aq (and hence ay) is large enough. 
Note that the estimates above are uniform in z varying in a small interval. It remains to 
show that 

Af 



29 



is uniformly integrable and bounded in z for z varying in a small interval. To this end we 
use that (f89|) is equal to 

EAf,y,L exp Gx, }p^,v,g.,l(q + ,■■■,«+ ''^ ^ 

j=i 

As in the proof of Theorem [5l we get 



P^y,G.,L(y) < C^' n Pkv,a,L(%), (90) 

where C is such that C > N/{N — k). By Corollary [25] in combination with Lemma [28] we 
have 

Pn,v,g.,l S L-e , 

where ||G.||oo := sup4g[_^^^] \Gt\- As ||G.||oo is sub-Gaussian, we get in combination with ([881) 
that (I86p is uniformly integrable w.r.t. P, provided that ag is large enough. It is clear that 
this bound is uniform in z G [0, e) for some small e > 0. 

To summarize, we have shown that (I85p converges to for (small) positive z, or in other 
terms, locally uniform convergence in a and t of VFAr,2(a, i) (for small positive z) as N ^ oo. 
We have also shown uniform boundedness of W^^ziO', t) for arbitrary A^, a, t and z £ (— oo, e) x 
M C C and as locally uniform convergence implies pointwise convergence, we get by Vitali's 
Convergence Theorem that the sequence (in N) of functions Wisi^zio-jt) converges to for 
z = —1 pointwise in a and t. To get locally uniform convergence in t and a for z = — 1, 
recall that by Arzela-Ascoli's Theorem, a sequence of continuous functions on a compact set 
has a uniformly converging subsequence if and only if the sequence is uniformly bounded 
and equicontinuous. Thus it remains to show that {Wiy^z{a,t))N' is equicontinuous in a and 
t (boundedness has already been shown). As the convergence of WN^zio-jt) is uniform in a,t 
for small positive z, Arzela-Ascoli's Theorem implies equicontinuity (in a,t) of {Wi\i^ziO',t))N 
for small positive z. To see that this implies equicontinuity (in a,t) of {WN,z{0',t))N also for 
z = —1, observe that a (real-valued) sequence of functions {fN)N on some compact K CM'^ 
is equicontinuous in x £ K if and only if for each sequence {xm)m C K, limm_i.oo Xm = 
X and each sequence {Nm)m C N we have Imim^oo fNmixm) — fNmi^) — 0- Using this 
characterisation, equicontinuity for z = —1 is easily seen by applying Vitali's Convergence 
Theorem to deduce limm^oo WN^-i{am,tm) = from limm_^oo WNmA'^rn,tm) = for smah 
positive z. This completes the proof of Theorem [3l 

To prove Theorem [H take g : — > M bounded and continuous. With the same arguments 
as above, we arrive in analogy to ([83]) - ([8^ at proving 



N 

^N,v,Lew{yZ^^j} / PN,V,G.,L(tl,---,tk)gitl,---,tk)dti...dtk 

N 

-EAr,i/,Lexp{V'Ga;.} / g{ti, . . . ,tk)tJ,iti) . . . fi{tk)dti . . .dtk 
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All the boundedness and integrability arguments above for En,v,lg^P {Ylf=i^Xj} can be 
used again. The convergence of f^t Pnvg {'t)s{'t)dt towards / g(t)ijL{ti) . . . /j{tk)dt is given by 
[Joh98[ Theorem 2.1]. Lemma [281 enables us to transfer Johansson's result to the correlation 
function p%yQ i- This finishes the proof of Theorem [H □ 

7. Appendix: Equilibrium Measures with External Fields 

In this appendix, we recall some results about equilibrium measures, mainly from the book 
by Saff and Totik |ST97t Section I.l]. The following can be found in [ST971 Section I.l]. 
Let A^^(S) denote the set of Borel probability measures on a set S. Define for S C C compact 
the logarithmic energy of ^ ^ as 

:= log\z - dfi{z)dfj,{t) (91) 



and the energy V of hy V := inf^g_A4i{s) Ht^)- It turns out that V is finite or oo and in 
the finite case there is a unique measure uj-£ which minimizes (I9ip . This measure ui-£ is called 
equilibrium measure of S and the quantity cap(S) := e~^ is called capacity of E. For an 
arbitrary Borel set S we define the capacity of S as 

cap(S) := sup{cap(i^') : K C T, compact}. 

Lemma 30. If T, = [—1,1], I > 0, then cap(S) = 1/2 and the equilibrium measure is the 
arcsine distribution with support [—1,1]: 

dujj:{t) = ^^=dt, te[-l,l]. 
vrv t — i 

loy; has mean and variance P/2. 

Proof See [ST971 Section I.l]. □ 

Definition 31. Let S C M be closed. Let Q : T, — > [0, oo] satisfy 

a) Q is lower semicontinuous, 

b) Sq := {t S S : Q{t) < oo} has positive capacity, 

c) if E is unbounded, then lim Q(t) — log \t\ = oo. 

If Q satisfies these properties, we call it external field on S and W = e~^ its corresponding 
weight function. 

Furthermore, define for ^ G A^^(S) the energy functional 

Iq{iA-= [Q{t)dfi{t)+ [ [log[s-t[-Ufi{s)dfiit). (92) 



Remark 32. In |ST97j the authors define the energy functional to be (in our notation) I2Q 
instead of Iq. It is more convenient for our purposes to use this definition. We note that 
under this change qualitative results from [ ST97] remain the same but quantitative results 
involving Q have to be changed by a factor 2 or 1/2, respectively. 

Iq{^j) might be 00, but the following theorem holds. The support of a measure ^ will be 
denoted as supp(/u). 



31 

Theorem 33. Let Q be an external field on S. 

a) There is a unique probability measure fig e with 

b) HQ has a compact support. 

c) Let Q be an external field on S such that Q = Q on a compact set K with supp(/iQ) C K 
and Q{t) = oo for t ^ K. Then fig = fiQ. 

Proof. Statements 1) and 2) can be found in |ST97t Theorem 1.1.3], 3) follows from [ST971 
Theorem 1.3.3] (also see the remark on page 48 in [ST97j ). □ 

fj.Q is called the equilibrium measure for Q. The next theorem summarizes properties of the 
logarithmic potential 



U^'{z) := j log\z-tf^ d^it). 



Theorem 34. 

a) Let Q and Q be external fields on S such that \Q — Q\ < e on S. Then for all z ^ C 

|C/^«(z) - U^Q{z)\ < 2e. 

b) Let K C be compact and {fin)n be a sequence in Ai^{K) converging weakly to a 
probability measure fi. Then for a.e. z G C (w.r.t. the Lebesgue measure on C) 

liminf C/^"(z) = C/^(z). 



71— >00 



c) // fi and V are two compactly supported probability measures and their logarithmic 
potentials and coincide almost everywhere on C, then = v. 

Proof. Statement 1. is contained in |ST97t Corollary 1.4.2], statement 2. is |ST97t Theorem 
1.6.9] and assertion 3. is [ST971 Corollary II.2.2]. □ 

We also need a characterization of the support of the equilibrium measure. 

Theorem 35. Let Q be an external field on S. 

a) For a compact set K of positive capacity define the functional 

Fq{K) := logcap(i^) -2 j Qdcox- 

For any compact K of positive capacity we have Fq(K) < i<Q(supp(^Q)). Furthermore, 
if K is compact and of positive capacity and such that Fq{K) = Fq(supp(/xq)), then 
supp(/xq) C K. 

b) If Q is convex, then supp(//Q) is an interval. 

c) If Q is even, then supp(//Q) is even. 

Proof. For statement 1. see [ST971 Theorem IV. 1.5], for statements 2. and 3. jST97t Theorem 
IV.1.10]. □ 

The last fact is about existence and properties of a continuous density of the equilibrium 
measure. 
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Theorem 36. 

a) Let Q he an external field on S. If Q is finite on supp(/iQ) and locally of class C^"^^ 
for some e > (which means that Q is continuously differentiahle and the derivative 
Q' is Holder continuous with parameter e ), then has a continuous density on the 
interior o/supp(/iQ). 

b) If Q has two Lipschitz derivatives and is strictly convex, then supp(^Q) =: [a,b] and 
the density of fiQ can be represented as 

^ = r(t)V(t-a)(6-i)l[,,,](t), (94) 

where r can be extended into an analytic function on a domain containing [a, b] and 
r{t) > for t G [a, b]. In particular, the density is positive on (a, b). 

Proof. Statement 1. is |ST971 Theorem IV.2.5], for assertion 2. see e.g. the appendix of the 
paper by McLaughhn and Miller |MM08| . □ 
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